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Abstract. - We generalize Schnakenberg's theory to nonsteady states, on networks representing 
currents between microstates of a nonequilibrium system, with the introduction of a new set of 
macroscopic observables which, for planar graphs, are related by a duality. We apply this duality 
to the linear regime, yielding a dual proposition for the minimum entropy production principle, and 
to discrete electromagnetism, finding that it exchanges fields with sources. We interpret duality 
as acting between system and environment, and discuss generalization to nonplanar graphs. 
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Introduction. — In a seminal paper [T], J. Schnaken- 
berg engaged in the definition of the fundamental macro- 
scopic observables of NonEquilibrium Statistical Mechan- 
ics, grossly conceived as a theory of the internal flows of 
a system. In accordance with modern trends in Quan- 
tum Field Theory and its outskirts, where Wilson loops 
and adiabatic phases play an ever more prominent role, he 
foresaw circuitations of certain "connection" variables as 
the constraints which prevent a system from relaxing to 
equilibrium. Born out of the study of biophysical systems 
[2,3, and recently finding growing applications to chem- 
ical reaction networks, molecular motors and transport 
phenomena @H9], his analysis has a deep geometrical and 
combinatorial content [10l[TT]. It is the backbone to the 
comprehension of Non-Equilibrium Steady States (NESSs) 
[IUJ[T2], to which the theory is restricted — so far. 

The aim of this letter is to go beyond NESSs, generaliz- 
ing Schnakenberg's construction to arbitrary states. The 
complete theory of nonequilibrium observables turns out 
to enjoy a duality which exchanges the concept of steadi- 
ness with that of detailed-balancing of the external forces. 

Duality comes in many flavours in physics. Among the 
first that one encounters: the duality between vectors — 
velocities — and linear forms — momenta; the Legendre 
transform which maps the lagrangian into the hamilto- 
nian, pivoting on the bilinear form Yj%QiPu the electro- 
magnetic duality, which is the archetypical physical coun- 
terpart of Hodge's geometrical theory of differential forms; 
and, less credited, the electro-technical duality between re- 
sistances and condensators, parallel and series reduction, 



voltage and current laws [13] ■ The one that we put for- 
ward descends from the latter, abstracting and general- 
izing it to nonlinear regimes, where Ohm's law does not 
necessarily hold; but it also resonates with each of the 
above. While the reference physical situation is that of a 
nonequilibrium system, we will cast our propositions in a 
very general form. In fact, they can be applied to any dis- 
cretized theory which has a couple of conjugate variables. 

Duality can only be realized on planar graphs. Al- 
though, nonequilibrium observables behave "as if" there 
always existed some dual graph. In a fascinatig work 
[13] . McKee attempts a generalization, finding a corre- 
spondence with logical duality between the universal and 
existential quantifiers (V and 3) under the involutive ac- 
tion of negation (->). In the prologue he comments that 
"some optimists see them [dualities] as mechanically dou- 
bling the number of results of a theory" . We enroll the 
troop, claiming that for every proposition that is true of 
steady states, there exists a dual proposition regarding 
detailed-balanced systems. One explicit example is the 
following dual minimum entropy production principle: 

"detailed balanced systems are those systems 
for which the rate of entropy production has 
the minimum value that is consistent with the 
fixed inflowing currents which prevent them from 
reaching a stationary state" . 

While steadiness is a property of the state of the system, 
detailed-balancing of the external forces is a property of 
the state of the environment: hence the letter's title. 
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Schnakenberg's main focus was on rate equations 

Pv = X i W vv'Pv' ~ w v , v p v ) = Yjiw, (1) 

v' v' 

where vertex v belongs to a finite state space V of the 
system, p is a probability density and w vv > are positive 
transition rates along edges e = v *- v' of a network, or 
graph, G. He identified the entropy production (EP) 

<t = Y, Wvv'Pv' In Wvv — = \Yj Jwa-vv' (2) 

v,v' W v ' v p v I v v , 

as a bilinear form of mesoscopic currents j vv > = -j v ' v and 
forces a vv > = -a v , v , and realized that at a NESS, that is, 
when Kirchhoff's law p = is satisfied, currents can be 
expressed as linear combinations of a certain number of 
macroscopic "mesh" currents J a . On the p = shell, EP 
comes down to a = J a A a (repeated indices are implicitly 
summed over), where the conjugate variables A a are seen 
to be circulations of the mesoscopic forces around suitable 
cycles of the graph, 



A a = In ■ 



w„ 



. w 



V3V2 w V2V1 



(3) 



A system whose steady state p ss makes all mesoscopic cur- 
rents and forces vanish, w vv 'p s J = w v > v p s v s , is said to satisfy 
detailed balance. Schnakenberg's choice of circuitations as 
the fundamental obscrvables, indicative of the noncquilib- 
rium nature of the system, is motivated by the well-known 
fact that they all vanish if and only if the steady state 
is detailed balanced (Kolmogorov's criterion) [12]. More- 
over, they do not depend on the system's macrostate p: 
they are external properties, which are conceptually more 
fitly awarded to the state of the environment. Hence, in 
the following, we will refer to detailed-balanced systems 
as those which satisfy Kolmogorov's criterion. 

Indeed, Schnakenberg's analysis can be extended to any 
graph whose edges bear a couple of antisymmetric con- 
jugate variables, one of which obeys Kirchhoff's Law at 
the nodes. Thence abandoning rate equation thermody- 
namics — but retaining the nomenclature, we complement 
Schnakenberg's definitions with a new set of conjugate 
macroscopic observables. 

Schnakenberg revisited. — The results are based on 
a decomposition theorem of the EP in cycles and flows (or 
cocycles) of the graph. To give a first hint, consider the 
3-level system depicted with straight lines in figUJi, 



J3,a 3 
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(4) 



By the Handshaking lemma Q£„ = 2£ e ), EP ([2]) can be 
recast as a = a%ji + a^ji + o<333- We reshuffle, add and 
subtract terms to obtain 




Fig. 1: (a.) A planar graph and its dual, with vertices de- 
picted as bullets, edges (arbitrarily oriented) by straight lines, 
dual vertices with boxes, dual edges with curved lines. Span- 
ning trees are solid, their complements are dashed. Shaded 
angles indicate clockwise orientations, (b.) The oriented cycle 
generated by ei, and its dual cocycle. (c. and d.) In straight 
lines, a foundamental set of cocyles, respectively generated by 
cochords el and e%,\ in curved lines, their dual cycles. 



where we defined, along with one Schnakenberg circulation 
A\ (figOi>), the macroscopic currents J| , flowing out of 
vertex w 3 (figflb), and Jg, flowing into vertex w 2 (figHtO- 
Since, by eq.([T]), J 2 = ~f>3 and J3 = f>2, it is conceptually 
appropriate to ascribe these observables to the state of 
the system. The vanishing of A\ provides balancing, the 
vanishing of J 2 ,J 3 defines steadiness. In graph-theoretical 
language, amd Jjj" are weighted cocycles, that is, edge 
sets whose removal disconnects the vertex set V into two 
noncommunicating components: they measure the total 
flow from one set towards the other. 

Generalizing, let G = (V,E,d) be a connected graph 
on \V\ vertices and \E\ edges. Edges carry an arbitrary 
orientation (a choice of tip and tail vertices), with -e des- 
ignating the inverse edge. With the exception of loops 
and multiple edges (which are nevertheless allowed), the 
topology of the graph is described by the incidence matrix 



+ 1, if «- v 
-1, if 5-v 
0, elsewhere 



(6) 



(ai + a 2 + a 3 ) ji + a 2 (32 a 3 (33 -ji), (5) 



We use a rather algebraic approach to graph theory 
15,16 , working with integer linear combinations of edges 
in the lattice £ = 1, E , upon which d acts as a boundary 
operator. It is a standard result that d induces an or- 
thogonal decomposition of £ = C @ C* into the cycle space 
C = ker(d) and the cocycle space C* = rowspace(d). The 
dimension of the cycle space is given by the cyclomatic 
number C = \E\-\V\ + 1, whence by the rank- nullity theo- 
rem the cocycle space has dimension |V) — 1. 

From a graphical point of view, cycles c are chains of 
oriented edges such that each vertex is the tip and the tail 
of an equal number of edges (possibly none) . It is simple 
if it is connected, has no crossings or overlapping edges. 
A simple cycle can exist in two opposite orientations. A 
simple cocycle c* is a collection of edges whose removal 
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disconnects the vertex set into two components; it might 
carry one of two possible orientations when all edges point 
from one of the two components, called the source set 
s(c*), towards the other (see fig|Tb,d). 

Of all possible integral basis of £ , we concentrate on 
fundamental sets, which are so built. Let T £ E be a 
spanning tree of the graph (i.e. a maximal subset of E 
containing no cycles); we call its edges e* the cochords. 
The remaining edges e a e E \ T are called chords. There 
are |V| - 1 cochords and C chords. When adding a chord 
to a spanning tree, a simple cycle c a is generated, which 
can be oriented accordingly with e a (see figOJi,b). The 
fundamental set of cycles {c a } so generated is a basis for 
C. Similarly, when removing a cochord e*, the spanning 
tree is disconnected into two components, which identify 
a simple cocycle c*, with orientation dictated by e* (see 
fig|T^,c,d). Again, the fundamental set of cocycles {c*} is 
a basis for C* . 

The crucial peculiarity of fundamental sets is that no 
chord is shared by two cycles, and no cochord is shared by 
two cocycles. Moreover, any of the sets {e Q , e* }, {c Q , e*}, 
{c",c*}, {e Q ,c*} forms a basis for £. 

We finally introduce: (I) mesoscopic currents j e as real 
edge variables, antisymmetric by edge inversion j_ e = -j e ; 
(II) antisymmetric conjugate forces a e ; (III) the EP 

a = EieOe = (j,a), (7) 

e 

where the r.h.s. is the euclidean scalar product on the 
edge set in shorthand; (IV) the macroscopic observables 

J a = (e a ,j), J£ = (c£,j), (8a) 

A a = (c a ,a), A; = (e;,a). (8b) 

In left-to-right, horizontal order: internal currents flow 
along fundamental chords, external currents are the total 
flow out of the source set of a cocycle, external forces are 
circuitations of forces along the fundamental cycles, inter- 
nal forces are exerted along edges of the spanning tree. 
In terms of the incidence matrix, KirchhofFs Law reads 
dj = 0, and the r.h.s. of eq.(fT]) reads p + dj = 0. 

Oriented overlappings between edge sets can be suc- 
cinctly expressed in terms of the scalar product: 

(c a ,e p )=6$, «, (9a) 

(c Q ,<)=0, (e;,e a )=0. (9b) 

At last we can formulate the core theorem, whose proof is 
given in the appendix. 

Theorem 1. The entropy production can be decomposed 
into a steady-state and a detailed-balanced term, a = a ss + 
o~db, given respectively by 

a ss = A a J a , a db = J^A;. (10) 

A detailed-balanced system (resp. steady state) has 
vanishing external macroscopic forces A a = (resp. cur- 
rents J* = 0), in which case o~ ss (resp. o~db) vanishes for all 
values of the internal currents (resp. forces). When both 
vanish we talk of equilibrium states. 



Duality. — A graph is planar if it can be drawn on the 
surface of a sphere with non-intersecting edges. Planar 
embeddings have faces / 6 F, i.e. open neighbours of the 
sphere which cannot be path-connected without crossing 
an edge. Their number = C + 1, including the "outer" 
face, is prescribed by Euler's formula. 

The dual graph G* = (V*,E* ,d*) has one vertex per 
face, V* = F, two dual vertices being connected by one 
dual edge *e per each boundary edge e that the corre- 
sponding faces share, so that E* = E. Pictorially, after 
puncturing and flattening the sphere, one will draw a ver- 
tex inside each face and a dual edge *e crossing e, then 
assign an orientation by clockwise rotating *e until it over- 
laps, tip and tail, with e (see curved lines and shadings in 
figfTJi). Crucial facts about duality are: 

(i) Up to a reorientation E -*■ —E, it is involutive; 

(ii) Different embeddings might yield non-isomorphic du- 
als (with different incidence relations); 

(iii) It maps a spanning tree T to the complement T* = 
E\T* of a spanning tree T* £ E* , in such a way that 
the fundamental sets generated by T* are the duals 
of the fundamental sets generated by T, according to 
the scheme (see fig|Tb,c,d) 

chords ** cochords, 
cycles <-»• cocycles. 

Duality can then be applied to the graphical structure 
of nonequilibrium observables. So, for example, the map 
a ** j leaves a invariant, but switches macroscopic observ- 
ables with those of the dual graph, mapping internal (resp. 
external) forces to internal (resp. external) currents: 

A a ±+ Ji\ A\ L ~ J a , <T SS <T db . (11) 

Since we ascribed A a to the state of the environment 
and J* to that of the system, it is fair to dub this 
system-environment duality. Steady states, for which the 
macroscopic external currents vanish, are dual to detailed- 
balanced systems, for which the macroscopic external 
forces vanish: the former are in fact properties of the sys- 
tem under given enviromental conditions, while the second 
are properties of the environment's influence on the sys- 
tem, indepently of the system's state. 

Out of the a ** j special case, we stress that duality 
is a graph-theoretical property: it tells how well-behaved 
observables should look like from the point of view of the 
environment and of the system, not which mesoscopic vari- 
ables should play the game. 

Planarity seems to be a major limitation to the general- 
ity of system/environment duality. We argue that this is 
not the case. Property (iii) listed above is independent of 
the particular embedding chosen. Indeed, generalizing the 
concept of a graph to that of an abstract matroid [TTIHB] , 
it turns out that matroids always have a well-defined dual 
which satisfies property (iii), even though dual matroids 
might not be visualizable as graphs. 
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Linear regime and minimum entropy produc- 
tion. — One major clue that led Schnakenberg to iden- 
tify chords and cycles as good thermodynamic observ- 
ables is the fact that, in the linear regime, Onsager's reci- 
procity relations arise. By "linear regime" it is meant 
that mesoscopic currents and forces satisfy Ohm's law 
a = £j + 0(j 2 ), where £ = diag{^i, . . . ,£\e\} is a "local" 
linear response matrix, connecting mesoscopic quantities 
edge by edge. Suppose that a system is initially at equi- 
librium, and it is perturbed to a nearby nonequilibrium 
steady state. Schnakenberg furnished the macroscopic lin- 
ear relation A a = L a ^Jp 1 with a symmetrical linear re- 
sponse matrix L. In our algebraic formalism, the deriva- 
tion is very strightforward: 

A a = (c a ,£j) = {c a M).h = L a0 Jp- (12) 

The linear response matrix is a weighted superposition 
of cycles. For rate-equation systems, this insight is com- 
plemented by Andrieux and Gaspard's proof of a Green- 
Kubo-type of formula for L pQ. Let us now linearly 
perturb an equilibrium state into a nonsteady, but still 
detailed-balanced configuration. Like steadiness is implied 
by Kirchhoff's law, balancing follows from the dual rela- 
tion d*a = 0, which is solved by a = A*c*. Then 

J* = (c^o) = (c&r'O = V?K (13) 

and the dual response matrix L* is a weighted superposi- 
tion of cocycles. Both matrices L and L* are symmetric, 
and under £ V 1 are dual one to the other. Similar ma- 
trices are of use in electrical circuit analysis [T3] and in 
the parametric formulas for Feynman diagrams (see [19] , 
§18.4; [16], §3). In this contest planar-graph duality has 
been related to duality between momentum and position 
representations [20]. Possibly the most interesting prop- 
erty of L and L* is that their determinants, which are 
always nonnull but for very trivial graphs, are indepen- 
dent of the fundamental sets chosen, obey the relation 
detL/ dctL* = det£, and are related to the 0-state Potts 
partition function [21] . Another important fact is that 

Theorem 2. When the equilibrium state is linearly per- 
turbed in an unconstrained manner (neither into a steady 
state nor into a detailed balanced configuration), entropy 
production can still be written as a block- diagonal bilinear 
form of the external observables 

a = (■L- 1 ) aft A a Af > + (L^V-M- ( 14 ) 

This result is far from trivial — for example, cross-terms 
will appear when expressing EP in terms of the internal 
observables! — and further supports the point of view that 
the external currents and forces are the true macroscopic 
nonequilibrium quantities which the observer controls. It 
will deserve further attention both for its technical sub- 
tlety and for its conceptual simplicity. 

One physically-motivated application of Schnakenberg's 
macroscopic observables in the linear regime was proposed 



by the author [22], who proved that Schnakenberg's affini- 
ties are the correct macroscopic constraints to be imposed 
to the minimum entropy production principle, which in 
one particularly suitable wording [33] asserts that 

"the steady state is that state in which the rate of 
entropy production has the minimum value con- 
sistent with the external constraints which pre- 
vent the system from reaching equilibrium" . 

The key step is to identify the external constraints with 
Schnakenberg's external forces. Varying a(j) = (j,£j), 
with given values of (c a ,£j) = A a held fixed thorugh La- 
grange multipliers, 

^{o-(j)-2\ a [(c a ,£j)-A a ]} = 0, (15) 

we find the extremal solution j = X a c a . To fix the value of 
the Lagrange multipliers, replace in the constraint equa- 
tion to obtain L a °Xp = A a . Since L is invertible, A Q turns 
out to be the fundamental internal current conjugate to 
the observed value of the affinity A a : hence we are in 
a steady state. One can then easily dualize the whole 
derivation to prove the dual proposition quoted in the in- 
troduction. 

Electromagnetism on a network. An important 
notion of duality against which to compare ours is the 
elecromagnetic (EM) duality. We refer here to C.Timm's 
work on rate equations [24 . 

Let's think of p as a charge density. In order to make 
the overall network neutral we introduce a supplementary 
vertex "oo", charged poo = ~T,vPv- All graph-theoretical 
notation will refer to this extended graph, which can be 
further made into a 2-dim. cell complex by introduc- 
ing a collection P 2 C of plaquettes. Choose a clock- 
wise/anticlockwise orientation for each plaquette p and 
define the boundary (curl) operator (see |25] for an intro- 
duction on discrete differential calculus) 

f +1, if e IO P, e tQ p 
(dx)P = l -1, if etOp, elQp . (16) 
[ 0, elsewhere 

Boundaries of plaquettes (columns of dx) are cycles, hence 
d(dx) = (the divergence of the curl vanishes). 

Introduce an electric field E e over edges and a magnetic 
field Bp over plaquettes. The electric field is required to 
satisfy Gauss's law <9E = p. Taking the time derivative, 
we have d(E + j) = 0, from which it follows that E + j is a 
linear combination of cycles, 

i = -E + f a c°E-E + axB, (17) 

where the r.h.s. imposes the Ampere-Maxwell Law. Since 
c" is a complete set of cycles, there exists an |P|xC matrix 
rj such that (<9x) p = rf a c a , so that SS a = if a B v . Impose 
Faraday's Law (<9x) T E + B = 0, and apply rj: 

(c a ,E) = -m a . (18) 
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It follows that any two combinations of plaquettes which 
share the boundary enclose a volume across whose bound- 
ary the magnetic flux is zero. Only C out of |P| magnetic 
field values are independent (Gauss's Law). 

As EP it is reasonable to elect the total energy flux 

a = (E, E) + £ B p B p = -(j, E) = (E, E) + S§ a S8 a (19) 
p 

where we applied Faraday's Law, transposed the curl op- 
erator, and used Ampere's Law to get the second identity 
(Integrated Poynting's Theorem). The third displays a 
simple dependence on the boundary values of the mag- 
netic field. Our theorem can now be applied: 

a = j a m a - j>;g; (20) 

where S* is the electric field along cochord e* . By eq. (|2"2"j) , 
is (minus) the time-derivative of the charge in s(c*). 
Hence under graph duality and j E one obtains 

§* J a , g§ a ++p> i + const. (21) 

The electric field is mapped to the source of the magnetic 
field and viceversa. Thus the example further supports 
the interpretation of duality as reversing the role of sys- 
tem and environment. Although, notice that the dynami- 
cal evolution is not respected: only Kirchhoff 's and Fara- 
day's "structure" equations are dual to each other. The 
lagrangian turns out not to be self-dual. This is an impor- 
tant difference between sys./env. and EM duality, which is 
dynamical. Moreover, the former is 2-dimensional, while 
the latter, restricted to sourceless cases or requiring mag- 
netic charges, involves the Hodge machinery in 3 dimen- 
sions. In ours, divergencelessness of the magnetic field is 
an essential feature rather than an obstruction to duality. 

Discussion and conclusions. — The latter applica- 
tion highlights that duality only works for kinematical 
states, viz. istantaneous snapshots of the system. So, 
for example, by "steady" we mean that Kirchhoff 's law is 
satisfied, not persistence in time. This is one important 
limitation that one will have to take care of when con- 
sidering, for example, markovian evolution: by no means 
do we claim that duality maps the a rate equation into 
a dual rate equation. As to the other important limita- 
tion, namely planarity, we already discussed how it can 
be formally overcome with matroids and conceptually re- 
garded as accidental. However, there is another way out, 
based on the fact that any graph can be embedded on an 
orientable closed surface of high-enough genus, and there 
dualized. While the number of cycles and cocycles is not 
affected, the number of faces, hence of dual cocycles and 
dual cycles, will change according to Euler's formula. It is 
simple to foresee that Theorem Q] will hold unchanged, but 
its interpretation will have to be accordingly modified, ac- 
counting for a number of global currents and of topological 
phases, such as those which were taken in consideration by 
Jiang and the Qians [10]. Duality will only hold locally. 



To conclude, let us linger on the 3-state example, in the 
attempt to provide a more intuitive grasp on the physics of 
duality. Label the graph with the energy levels of an open 
system, which can emit and absorb enviromental quanta 
of energy. The onset of a NESS might be due to the in- 
teraction with two thermal baths 2(i 27 . whose inverse 
temperatures Pa and Pb label the states of the dual sys- 
tem, with Pa > Pb- Suppose that transitions 2 and 3 are 
exclusively due to the interaction with B, while transi- 
tion 1 is exclusively due to the interaction with A, as is 
pictorially depicted in figG^. The ratio of emission and 
absorption rates is given by w ei /w- ei = exp@A(v2 _ W3), 
and similarly for the others, yielding as macroscopic affin- 
ity A\ = (/3a - Pb)(v2 ~ V3). In a nonequilibrium steady 
state, with current j± = j 2 = J3 = J\, one quantum of en- 
ergy v 2 - v 3 is emitted towards reservoir A every |Ji| _1 
seconds, per each couple of quanta of energy v 2 - v\ and 
v\ - V3 that are absorbed from reservoir B. It takes shape 
a picture where to a steady state there corresponds a non- 
steady flow of energy from the hotter to the colder bath: 

nonequilibrium steady sys. -*■ nonsteady env. 

Whilst purely speculative, this interpretation is consistent 
with the physical intuition that NESSs are determined by 
a transient enviromental behaviour [28] . 

Viceversa, a detailed-balanced flow arises when there is 
no temperature gradient, Pa = Pb, in which case we only 
resolve one reservoir. At equilibrium, because of steadi- 
ness and detailed balancing, as many quanta of given en- 
ergy are emitted and absorbed. However, fluxes within 
the system determine a non-null flow of currents in the 
one bath, which is necessarily steady being a 1-state sys- 
tem. Hence the system's state plays the role of external 
force igniting internal fluxes within the environment: 

steady env. -> detailed balanced nonsteady sys. 

This is nothing but the logical negation of the above 
proposition, hence its dual under transposition of the ma- 
terial implication symbol (->•). 

Despite of its simplicity, the example is rather clumsy 
and only vaguely illustrative: system and environment do 
not play mirror roles, for which reason we were not able 
to draw the inverse implications. However, the qualitative 
principle seems to be robust. It is quite remarkable that 
graph duality finds a similar interpretation also in me- 
chanical engineering [29] , where the statics of a structure 
and its first order kinematics are related to dual proper- 
ties of its design. Thus there seems to be a vast variety 
of systems to which duality might apply: it is the au- 
thor's opinion that the development of a statistical model 
which displays duality between enviromental and internal 
degrees of freedom would be a major advance. 

* * * 

The author warmly thanks A. Maritan and D. Andrieux 
for discussion, M. Esposito and M. Dalmonte for helping 
out with the first drafts, and F. Sbravati for support. 
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Fig. 2: (a.) Transitions between states due to absorption and 
emission from two reservoirs, (b.) Steady state heat trans- 
fer between reservoirs, (c.) One reservoir with internal fluxes 
ignited by the system's nonsteady configuration. 

Appendix. 

Proof of Theorem [7J The strategy is to find the general 
solution to the continuity equation with sources p+dj = 0. 
Here p e R v is the current injected at the vertices, and it 
is constrained to satisfy Y*vtV Pv = 0. Since any |V| — 1 rows 
of d span the cocycle space, p v is expressible as a linear 
combination of a fundamental set of external currents, and 
viceversa. One can easily show that 

-P M = J£ = - E P* (22) 

Physically: the flow out of a source set is equal to (minus) 
the sum of the injected currents within the set. 

The general solution can be found as a particular solu- 
tion plus the general solution of the homogeneous equa- 
tion associated to it. Solving dj = yields a superposition 
of cycles Z a A Q c Q . As to the particular solution, since 
{c",e*} is a basis for £, we can tune the cycle currents 
so as to make currents along chords vanish. We then only 
need to specify the particular solution along cochords, ob- 
taining 

j = X a c a + A{fe*. (23) 

Inserting (|2U)) into the definitions (|8ap . and using the or- 
thonormality relations one identifies J a = X a , and 
J* = A*. Further insertion into ([7]) yields our thesis. □ 

Proof of Theorem^ Consider eq.(|23j). with X a = J a and 
A* = j!f, and replace in the bilinear form a = 

a = L aP J a Jp + M^JfJ" + 2Hp a J^ (24) 

where we defined Af" = (e£,&£) and H% = (c a ,£eJA. 
Moreover, it's simple to calculate A a = L a @ Jp + H^J*. 
Completing the square: 

a = L^A a Af + (M M „ - HZL^HftJW. (25) 

Since A a and Jf are independent, setting all affinities to 
zero yields the entropy production for detailed balanced 
systems, which after the previous theorem and eq.(|13j) is 



easily seen to be a = (L^ 1 )^ J* J". Since the latter is a 
nondegenerate bilinear form, we can identifty the matrix 
between parenthesis with L^ 1 . □ 
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